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1 . INTRODUCTION 


The  search  for  models  useful  in  describing  and  predicting  the 
break-up  of  shaped-charge  jets  has  proceeded  along  several  lines. 

The  problem  of  mass  flow  in  jets  can  be  treated  in  a continuum  mechanics 
fashion  by  formulating  a set  of  difference  equations  which  can  be  solved 
in  a reasonable  time  by  a large  computer.  When  spatially  periodic 
perturbations  are  introduced  into  the  boundaries  of  the  jet,  they  are 
found  to  grow  if  they  are  of  the  right  size.*  Another  approach 'is  to 
make  experimental  observations  and  attempt  to  use  similarity  methods2 
and  empirical  curve  fitting  to  extrapolate  from  the  measured  to  the 
unknown.  A third  approach,  and  the  one  to  be  followed  in  this  study, 
is  to  use  basic  mechanics  and  simple  probability  theory  to  treat  the 
jet  as  a collection  of  stretching  links.  All  of  these  approaches  are 
tied  to  material  properties  which  are  reasonable  in  concept  but  dif- 
ficult to  measure  during  the  formation  and  break-up  of  the  jet 

It  has  been  observed  that  high  speed  jets,  created  by  the  detonation 
of  shaped  explosive  charges  with  conical  copper  liners,  tend  to  break 
up  into  small  particles  of  various  sizes  after  about  10”^  seconds. 
Continuous  jets  with  all  parts  traveling  along  the  same  straight  line 
are  known  to  have  better  penetration  capabilities  than  jets  which  are 
broken  or  dispersed,  and  there  is  considerable  interest  in  discovering 
and  controlling  those  physical  parameters  which  influence  the  pattern 
of  jet  break-up.  The  purpose  of  this  report  is  to  combine  a number  of 
observations  with  some  simple  theory  and  to  arrive  at  a model  which 
contains  the  pertinent  physical  parameters  and  predicts  the  pattern 
of  the  break-up  of  shaped  charge  jets. 


2.  SIMPLE  DESCRIPTION  OF  A JET 


By  examining  consecutive  x-ray  flash  photographs  of  a jet,  it  is 
possible  to  determine  its  mass  distribution  in  space  and  time  and  to 
describe  it  rather  simply  at  a given  time  as  shown  in  Table  I,  which 
will  be  the  basis  for  the  following  analysis.  Jet  properties  will  of 
course  vary  with  the  design  of  the  shaped  charge.  The  values  listed 
in  Table  I were  obtained  by  a crude  averaging  of  data  from  firings^ 
of  precision  3.3  inch  (8.4  cm)  shaped  charges  with  copper  liners, 
and  are  thought  to  be  representative.  Because  copper  is  the  metal 


P.C.  Chou  and  J.  Carleone,  Stability  of  Shaped  Charge  Jets,  Journal 
of  Applied  Physics,  Vol.  48,  No.  10,  pg.  4187,  October  1977 

“W.  E.  Baker,  P.  Wcstine  and  F.  T.  Dodge,  Similarity  Methods  in 
Engineering  Dynamics,  Hayden  Book  Co.,  Inc.,  Rochelle  Park,  NJ , 
pg.lhf),  1973 

’.I.  Simon  ami  R.  Ui  Persia,  Shaped  Charge  Warhead  Performance,  Transac- 
tion', of  4th  Symposium  on  Warhead  Research,  US  Naval  Ordnance  Testing 
'mat  ion,  China  Lake,  CA,  Sept.  1963 
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most  often  used  in  liners  for  shaped  charges,  some  properties  of  copper 
are  listed  in  Table  II. 

-4 

Table  I.  Nominal  Properties  of  a Copper  Jet  10  Seconds 
After  Shaped  Charge  Initiation 


Velocity  of  Tip 
Length  (V>2.5  km/sec) 
Diameter 
Break-up  Time 
Number  of  Particles 
Mass  in  Jet 


7.5  km/sec 
50  cm 
3 mm 

100  u/sec 
50 

32  gms 


Table  II.  Nominal  Properties  of  Copper 


Density 

■5 

8.93  gm/cm 

Sound  Velocity 

3750  m/sec 

Thermal  Conductivity 

3.94  w/cm°C 

Heat  of  Fusion 

207  j/gm 

Heat  of  Vaporization 

4730  j/gm 

Specific  Heat 

0.38  j/gm°C 

Atomic  Weight 

63.54 /mol 

Melting  Point 

1083°C 

Boiling  Point 

2595°C 

Strength 

2X108  newton 
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For  purposes  of  discussion  we  will  assume,  as  depicted  in  Figure 
1,  that  the  shaped  charge  produces  at  the  virtual  origin  a copper  jet 
which  expands  uniformly  in  time  so  that  it  can  be  viewed  as  a stretching 
right  circular  cylinder  of  constant  mass  and  density.  The  length  L 
of  the  jet  at  time  t is  given  by  the  equation, 

L = Ugt  * 5.0  X 103  t (2.1) 

where  Ug  is  a stretching  speed  of  5,0  km/sec,  the  difference  in  speed 
between  the  tip  and  the  tail  of  the  jet.  If  the  jet  has  a diameter  D 
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Figure  1.  A depiction  of  the  jet  resulting  from  a copper-lined 

shaped  charge.  If  we  assume  that  the  jet  started  from 
the  virtual  origin  at  time  zero,  it  is  100  microseconds 
old.  The  speed  of  stretching  l),  is  5,0  km/sec  and  the 
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velocity  gradient,  or  strain  rate,  »s/L,  is  10  sec 
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of  3 mm  when  it  is  50  cm  long  at  100  p/sec  and  if  copper  has  a density 

3 

of  8.93  gm/cm  , the  mass  of  the  jet  is  given  by  the  equation, 

M = itD2Lp/4  = 31.56  gm.  (2.2)  ' | 

The  motion  of  the  jet  can  be  described  in  terms  of  the  motion  of  the 

center  of  mass  and  the  motion  relative  to  the  center  of  mass  by  the  ' 

equation, 

‘ i 

U(X)  = Uc  + (X/L)  Us;  - L/2<X<L/2  (2.3) 

where  U(X)  is  the  speed  of  a mass  point  located  by  the  coordinate  X 
measured  from  the  center  of  mass,  Uc  is  the  speed  of  the  center  of 

mass,  L is  the  length  of  the  jet  at  time  t,  and  U£  is  the  difference  j 

in  speed  between  the  jet  tip  and  tail.  | 

Combining  Equations  2.1  and  2.3  yields  the  result,  ! 

: I 

I 

U(X,t)  = Uc  + (X/t),  -Us/2<X/t<yjs/2,  (2.4)  J 

where,  for  our  particular  jet,  both  Uc  and  U are  5 km/ sec. 


3.  ENERGY  DISTRIBUTION  IN  THE  JET 


It  is  evident  that  Equation  2-4  can  hold  only  if  the  stretching  of 
the  jet  requires  no  expenditure  f energy.  The  energy  of  the  jet  may 
be  divided  into  two  parts,  the  one.’gy  associated  with  the  motion  of 
the  center  of  mass  and  the  energy  associated  with  motion  relative  to 
the  center  of  mass,  and  may  be  written  as 


E 


MU2/ 2 


d/2) 


r 


UU  (X)dX 


X=-L/2 


(3.1) 


where  p is  the  mass  per  unit  length  in  the  jet  and  'u(X)  is  defined 
from  Equation  2.3  by  the  equation, 


u (X)  = U(X)  - Uc  = (X/L)US. 


After  integration  between  limits  Equation  3.1 


uocume 


(3.2) 


I 


i 


i 
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(3.3) 


E = MU2/ 2 + MUg/24 

where  we  have  replaced  p by  M/L.  Since  Ug  and  U£  are  both  equal  to 

5.0  km/sec  for  our  situation,  we  see  that  only  about  1/13  or  8%  of 
the  energy  is  associated  with  the  stretching  motion,  and  it  is  this 
energy  which  is  depleted  as  work  is  done  in  stretching. 

Suppose  that  a uniform  jet  is  divided  into  n equal  segments  with 
the  velocity  distribution  of  the  Kth  segment  relative  to  the  laboratory 
system  being  given  by  the  equation, 


U(X)  * Uc  + (X/L)U  , A<X<B,  (3.4) 

where 

A i (-L/2)  + (K-l)  (L/n),  (3.5) 

and 

B = (-L/2)  + K(L/n) . (3.6) 

With  these  definitions,  the  location  and  speed  of  the  center  of  mass 
of  the  Kth  particle  are  given  by  the  equations, 

XK  = -L/2  + (K-Js)  (L/n)  = (A+B)/2,  (3.7) 

and 

UK  = Uc  + [-*5  + (K  -Jj)/n]  Us  - U(XK).  (3.8) 


The  velocity  distribution  in  each  segment  relative  to  its  center  of 
mass  can  be  written  as 


U" (X VL)  3 (X'/L)  Us,  -L/2n<XXL/2n 

(3.9) 

so  that  the  relative  kinetic  energy  can  be 

written  as 

X=L/2n. 

E'  = / (M/2L)  (X'/L)2  Us2  dX' 

(3.10) 

*/X=-L/2n 

L/2n 

3 (M/2L)  (Us/L)2[X'3/3] 

-L/2n 

(3.11) 

«=  MUg2/24  n3 

(3.12) 
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From  Equation  3.8,  the  kinet’*'  energy  of  the  Kth  particle,  less  its 
internal  rnergy  relative  to  its  center  of  mass,  can  be  written  as 


Fk  = (M/2n)  UjJ 


(3.13) 


- (u/2n)  (Uc  * [-h  + (k-y/n)U 


(3.14) 


Jsing  Equations  3.12  and  3.13  the  total  kinetic  energy  of  the  jet 
can  be  written  as 


K=n  K=n 

E » £ E,.  + £ E' 

K = 1 K K=1  K 


(3.15) 


K»n  7 K=n  - - 

£ (M/2n)  (U  + [-'-4  + (X-Jj)/n]U  CT  + £ MUV24nJ  (3.16) 

K=1  K-l 


Using  the  identities. 


K=n 

£ a = na,  a 2 any  constant, 
K=1 


(3.17) 


K~n 

= (*i)n(n+l) , 
K=1 


(3.18) 


K=n  - 

£ r = (1/6) (n) (n+1) (2n+l) , 
K-l 


(3.19) 


Equation  3.16  can  be  simplified  to  the  form, 


E = MU^/2  + (M/24)U^i.(i/n2)]  +(M/24)U*/n2 


= MU^/2  + (M/24)U^  . 


(3.20) 

(3.21) 


12 


Since  n segments  contribute  to  the  last  term  of  Equation  3.21,  tne 
contribution  of  each  segment  to  the  terra  is 

En  = (l/n)(M/24)(U2/n2)  = (M/24)U2/n\  (3.22) 

in  agreement  with  Equation  3.12.  Equation  3.20  states  that  the  n seg- 
ments contribute  to  the  total  energy  by  three  terms.  The  first  term  is 
associated  with  the  motion  of  the  mass  M at  speed  Uc.  The  second  term 

arises  from  the  motions  of  the  centers  of  mass  of  the  several  segments 
relative  to  the  center  of  mass  of  the  jet.  The  last  term  arises  from 
the  internal  motion  of  each  segment  relative  to  its  center  of  mass, 

The  total  energy  of  the  jet  is  of  course  independent  of  the  manner 
of  its  division,  as  illustrated  by  Equations  3.20  and  3.21.  These 
equations  can  bo  written  in  terms  of  energy  per  unit  mass  as 

E/M  = U2/2  + (l/24)ll2[l-(l/n2)  ] + (l/24)U2/n2  (3.23) 

and 

0/M  = U2/2  + (1/24)U2.  (3.24) 

The  above  calculations  illustrate  that,  with  a given  stretching 
speed,  the  energy  available  to  stretch  an  isolated  segment  of  jet, 
after  a partial  break-up  has  occurred,  varies  strongly  with  length. 

It  should  be  stated  further  that,  if  a jet  in  the  process  of 
stretching  suddenly  breaks  into  n pieces,  the  total  kinetic  energy  is 
still  the  same,  but  the  energy  available  per  unit  mass  to  stretch  the 

2 2 2 

jet  (or  its  several  pieces)  has  fallen  from  Ug/24  to  Us/24n  , as 
indicated  by  Equations  3.23  and  3.24. 


4.  RELATIVE  ENERGY  CONSUMED  IN  STRETCHING 

The  element  of  energy  dE  consumed  in  stretching  a rod  an  elemental 
distance  dX  is  given  by  the  equation. 


dE  ■ F dX,  (4.1) 

where  F is  the  applied  force.  By  the  definition  of  stress  S as  force 
per  unit  area  and  strain  de  as  elongation  per  unit  length,  Equation 
4.1  can  be  written  as 
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A X (dX/X)  = V S de, 
and  V - AX. 

. For  purposes  of  illustration,  suppose  that  the  jet  from  a 
shaped  charge  is  assumed  to  start  from  a right  circular  cylinder  of  some 
initial  length  Lq  and  that  it  continues  to  stretch  uniformly  until 

it  breaks  up  at  some  length  From  the  definition  of  de  Equation  4.3 

can  be  written  as 


(4.2) 

(4.3) 


dE  = (F/A) 
where  S = F/A,  da  e dX/X 
Then,  E/V  = J^Sde. 


E/V  = 


S (dX/X)  = S In  (Lf/L  ) 


(4.4) 


where  S has  been  assumed  to  be  a constant  parameter  during  the  stretching 
process . 

Suppose  that  a jet  exists  initially  as  a right  circular  cylinder 
with  a length  to  diameter  ratio  (Lq/Do)  of  unity.  The  volume  V of  the 

jet  can  be  written  as 


V = M/p  = it  0^  Lq/4  = ir  L3/4,  (4.5) 

so  that,  using  values  for  the  mass  M and  density  p from  Tables  1 and  II, 
the  initial  length  Lq  is  calculated  as 

Lq  = (4M/ifp)1/3  = 1.67  cm.  (4.6) 

Using  Equations  3.24  and  4.4  and  assuming  that  all  the  internal  kinetic 
energy  of  the  jet  is  used  in  stretching  and  that  there  is  no  jet 
break-up, 

Ug/24  * (S/p)  In  (Lf/LQ)  (*.7) 

and 

uf  = LoEXP[pU^/243].  (4.8) 
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= 1.64  EXP [8 . 93  x 103  • 50002/(24  • 2x  108) ] (4.9) 


2.64  x 1020  cm. 


(4.10) 


This  result,  though  perhaps  unrealistic,  illustrates  that  the  internal 
mechanical  energy  of  a jet  is  sufficient  to  stretch  it  until  it  breaks. 
If  we  had  assumed  a break  up  length  of  SO  cm,  the  work  done,  from 
Equation  4.4,  would  have  been 

E/pV  = (2  x 108/8.93  x 103)  In  (50/1.67)  (4.11) 

= 76.1  j/gm. 

This  is  only  about  7%  of  the  available  internal  kinetic  energy, 
given  from  the  last  term  of  Equation  3.24  as 

E/M  = U2/24  = 1042  j/gm.  (4.12) 

From  Table  II  the  specific  heat  of  copper  is  0.38  j/gm°C  and  the 
heat  of  fncion  is  207  j/gm,  sc  that  the  work  calculated  in  Equation  4.]] 
is  rather  modest  compared  to  the  energy  required  to  heat  copper  from 
500°C,  for  example,  and  melt  it.  Since  the  stress  S,  assumed  for  the 
foregoing  calculations  to  be  that  of  cold  copper,  is  expected  tc 
decrease  markedly  at  copper  approaches  the  melting  point,  there  is 
little  likelihood  that  work  done  on  a stretching  jet  will  cause  it  to 
melt  if  it  was  a solid  at  some  time  during  the  stretching  process. 

It  will  now  be  assumed,  in  agreement  with  observations,  that  the 
stretching  jet  will  break  rather  than  stretch  indefinitely.  In 
physical  terms  it  can  be  said  that  there  is  apparently  some  threshold 
of  work  by  stretching  which  the  jet  can  endure  before  it  breaks  up 
into  segments.  Each  segment  has  a velocity  gradient  U^/L  and  will 

continue  to  stretch  until  it  breaks  further  or  belief  waves  propagate 
from  the  free  ends  and  the  velocity  gradient  is  removed.  The  energy 
associated  with  the  velocity  gradient  persists  until  it  is  expended 
by  work  done  in  stretching.  The  energy  available  to  stretch  a segment 
of  length  l is  given  from  Equatioh  3.22  as 

E = (M/24)U2/n3  = (1/24)  (MJl/L)  (l^/L)2*2  (4.13) 

where  n has  been  replaced  by  L/i  and  factors  have  been  associated 
to  define  the  mass  Mfc/L  of  the  segment  and  the  velocity  gradient  Us/L. 
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The  internal  energy  per  unit  mass  in  the  segment  is  then 


K/M(d/L)  « (1/24)(Us/L)2«,2  • (4.14) 

Suppose,  for  example,  that  the  first  breaks  occur  at  some  energy 
threshold  ( E/M) Q and  that  subsequent  breaks  occur  if  the  segment 

has  enough  relative  energy  to  push  the  work  done  on  the  segment 
beyond  some  upper  limit  (E/M)^. 

Equating  the  additional  work  needed  for  further  break-up  to  the  energy 
available,  as  given  by  Equation  4.14,  we  obtain 

(1/24)(U  /L)V  = (E/M)  - (E/MV  s AE.  (4.15) 

b m u 

If  we  argue  that  the  segment  will  not  break  further  if 

(1/24)  (U  /L)V  < AE,  (4.16) 

then  the  equation, 

l < [24  (AE)/(Us/L)231/2  (4.17) 

defines  the  longest  segment  that  can  exist  without  further  break-up. 

By  means  of  Equation  4.4  and  4.15 


AE  = (S/p)  In  (Lra/l0),  (4.18) 

where  Lq  is  the  length  of  the  segment  when  break-up  starts  and  Em  is 

its  length  when  it  has  used  all  of  its  internal  kinetic  energy  in 
stretching,  Equation  4.17  can  be  written  with  Equation  4.18  to  yield 

« < [24  (S/p)  In  (Lm/L0)}l/,2/(Us/L)  (4.i9) 

Assuming  a 5%  stretch  (Lffl/Lo  = 1.05)  and  nominal  values  of  strength 
Sj,  density  p,  and  velocity  gradient  Ug/L  from  Tables  I and  II  and 
Figure  1,  we  find  that 

l < 1.62  cm  (4.20) 

Data,  obtained  by  Jameson  and  exhibited  in  Figure  7,  show  an 
estimated  maximum  length  of  1.98  cm. 
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One  way  of  looking  at  the  energy  or  time  interval  associated  with 
jet  break-up  is  to  balance  the  energy  available  for  stretching  a jet 
particle  against  the  work  done  in  stretching.  For  example,  a particle 
of  length  i,  velocity  gradient  Ug/L,  and  mass  M will  have  available  for 

stretching  at  most  the  energy 


E = (M/L)(Us/L)2(t3/24)  (4.21) 


as  given  by  Equation  4.13.  The  energy  per  unit  volume  is  therefore 

Ep/[M(S,/L)]  = pA2(Us/L)2/24  (4.221 

If  S is  the  strength  of  the  material  and  the  elongation  Ae  in  time  At 
is  given  by  the  equation 


Ae  - AA/i  = Us(£/L)  At/A,  (4.23) 

the  work  done  per  unit  volume  by  stretching  in  time  At  is  given  from 
Equation  4.2  as 


dE/V  « SAe  = SUs(i./L)  At/A.  (4.24) 

If  we  equate  the  enervv  available  to  the  work  done  in  time  At  (Equations 
4.22  and  4.24) 

pU2  (A/L)2/24  = SUs  At/L  (4.25) 


and 

At  = pUgL(A/L)2/24  S = p(Us/L)A2/24  S.  (4.26) 

If  we  use  the  value  of  A f~om  Equation  ,4. 20  we  find  that  (4.27) 

At  = 4.88  x 10"^  sec . 

From  purely  kinematic  considerations  using  Equation  2.1 

At  - t (AL/L)  = 5 X 10“6  sec.  (4.28) 


17 


The  time  internal  for  At  required  for  a relief  wave  to  travel  across  a 
1.62  cm  segment  is  given  by 

At  * 1/ a * 1.67  x 10"2/37S0  - 4.45  x 10-6  sec.  C4.29) 

As  will  be  shown  later,  this  transit  time  for  a relief  wave  is  intimately 
related  to  particle  size  distribution. 

The  point  of  the  preceding  discussion  is  that  maximum  particle 
length  (Equation  4.20)  and  relief  wave  speed  (Equation  4.29)  can  be 
consistently  related  by  a simple  stretching  model,  which  includes  work 
done  in  stretching. 

An  interesting  consequence  of  Equation  4.17  is  that,  if  the  amount 
of  work  permissible  by  stretching  becomes  small  in  the  time  interval 
between  when  break-up  starts  and  when  it  ends,  the  jet  will  break  into 
very  small  pieces.  The  theory  to  support  this  initial  conclusion  will 
be  developed  later. 

It  can  be  argued  that  the  assumptions  leading  to  Equation  4.20  and 
its  fair  agreement  with  experiment  are  ad  hoc  and  fortuitous.  It  will 
be  shown  later  that  the  same  assumptions  lead  to  a reasonable  distri- 
bution in  size  for  the  fragments.  Because  of  the  square  root  and 
logarithmic  dependence  of  i on  Lm/Lo,  the  result  is  not  particularly 

sensitive  to  fairly  large  variations  of  L /L  , For  example  if  L /L 

m o r m o 

had  been  10,  the  result  would  have  been  11.0  cm.  This  variation  in 
length  from  1.62  to  il  cm  represents  a factor  of  200  in  the  amounts  of 
stretch  assumed. 

The  velocity  gradient  (or  strain  rate)  Ug/L  is  probably  the  most 

easily  measured  of  the  quantities  in  Equation  4.19.  The  density  p is 
not  expected  to  vary  greatly.  The  strength  S,  on  the  other  hand,  is 
somewhat  of  a mystery  because  of  the  apparent  tendency  of  copper  to 
keep  its  strength  better  as  the  strain  rate  is  increased,  as  shown  in 
Figure  2. 

The  preceding  discussion  can  be  summarized  as  follows.  The  jet 
will  either  stretch  foreVet  (Equation  4.10)  or  it  will  break  up.  If 
it  breaks  up  and  its  time  of  break-up  is  determined  by  the  amount  of 
work  done  in  stretching  it,  there  will  be  some  threshold  (of  stretch, 
time,  or  work)  at  which  breakage  begins. 
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If  breakage  occurs,  each  segment  has  a velocity  distribution  and 
therefore  a kinetic  energy  relative  to  its  center  of  mass,  and  this 
kinetic  energy  must  be  consumed  in  further  stretching  of  the  segment. 
If  thcro  aro  further  potentially  weak  points  in  the  segment,  like  the 
ones  which  produced  the  segment  in  the  first  place,  the  segment  will 
subdivide  further  as  more  stretching  is  done.  This  process  will 
discontinue  rather  abruptly  because  of  the  sensitivity  of  the  energy 
available  for  stretching  to  the  segment  length  (Equation  4.14).  This 
process  can  be  modeled  simply  (Equation  4.19)  to  yield  an  estimate  of 
the  upper  limit  for  the  length  of  the  segments  following  break-up. 


5.  ELASTIC  AND  INELASTIC  BEHAVIOR  OF  A LONG  RIGHT  CIRCULAR  CYLINDER 
WITH  A LONGITUDINAL  VELOCITY  DISTRIBUTION 

In  order  to  better  understand  the  behavior  of  a jet  particle  with 
a velocity  gradient  and  two  free  ends,  it  is  useful  to  examine  the 
elastic  response  of  a long  right  circular  cylinder  to  a velocity  gradient 
as  an  initial  condition.  We  can  write  the  longitudinal  displacement 
y(X,t)  from  equilibrium  at  time  t of  a point  located  at  a distance  X 
from  the  center  of  mass  in  terms  of  the  Fourier  series, 


n=°° 

y(X,t)  = £ b sin  KnTrat/L)  + c ] sin|mrX/L)  + 5 ] (5.1) 

n n n 


where  L is  the  length  of  the  cylinder,  a is  the  speed  of  sound  in  the 
cylinder  and  bn,  en>  and  6^  are  constants  determined  by  initial  and 

boundary" conditions.  We  will  assume  that,  at  zero  time,  we  have 


\ 

y (X,0)  « 0,  -L/2CX<L/2, 

(5.2) 

and 

• i 

y(X,0) 

= UgX/L,  = RX,  -L/2CX<L/2, 

(5.3) 

i 

f 

i 

i.e.  that  the  rod  is 

initially  unstretched  but  has  ,a 

uniform  velocity 

j 

] 

gradient.  The  dot  over  the  y indicates  partial  differentiation  with 
respect  to  time,  Us  is  the  stretching  velocity  defined  oarlier,  and 

R ^ Us/L  (5.4) 

is  the  strain  rate. 
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The  compatibility  of  Equations  5.1  and  5.2  requires  that 

en  - 0.  (5.5) 

Differentiating  Equation  5.1  with  respect  to  time  t and  using  Equation 
5.5  gives  the  result, 

n“us 

y(X,t)  = 53  (nira/L)bncos  (nitac/L)sin[  (mtX/L)+6n] , (5.6) 

n=0 

so  that,  from  Equations  5.3  and  5.6, 


n="° 

y(X,0)  = 53  (nira/L)b  sin[(nirX/L)+6  ] = RX  (5.7) 

n=0  n 

The  Fourier  expansion  for  X in  the  range  -L/2<X<L/2  can  be  written'* 
in  the  form, 

]}SOO 

x = (4L/7T2)  £ [(-l)n'1/(2n-l)2]sin[(2n-l)xX/L]  (5.8) 
n=l 

Comparison  of  Equations  5.7  and  5.8  shows  that  it  is  possible  to  write 
the  coefficients  of  the  series  in  the  form 

6n  = 0 (5.9) 

1 

bn  ■ (4L2R/ir3a)  (-l)n-1/(2n-l)3  (5.10) 

Using  Equations  5.5,  5.9,  and  5,10,  Equation  1 can  be  written  as 

9 x n=“  „ l * , 

y (X, t)  = [4LTR/ir  a]  T [(-1)  /(2n-l)  ]sin[ (2n-l)irat/L]sin 

n=l 

[(2n-l)irX/L]  (5.11) 


“"Handbook  of  Tables  for  Mathematics  4th  Edition.  The  Chemical  Rubber 
Cor.  18901  Cranwood  Parkway,  Cleveland,  Ohio  44123,  pg.  611  formula  7, 
1964, 
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The  velocity  and  acceleration  are  computed  from  Equation  5.11  as 

n — °° 

y(x.t)  = C4LR/7T2)  £ [(-i)n+1/(2n-l)2]cos[(2n-l)rat/L]sin 

n=l 

[(2n-l)nX/L]  (5.12) 

and 

n=uo 

y(X,t)  = (4Ra/rr)  £ [ (-1)°/ (2n- 1) ]sin[ (2n- 1) irat/L]sin[ (2n-l) ttX/L]  (5.13) 
n=l 


Since  each  element  of  the  cylinder  is  accelerated  according  to 
Newton's  second  law,  we  can  write  for  the  mass  between  X and  X+dX  the 
equation,  dT  = pydX,  where  T is  the  tension  in  the  cylinder  and  u is  the 
mass  per  unit  length,  so  that 


(5.14) 


Using  Equat'on  5.13  and  the  fact  that  the  tension  at  the  end  of  the 
cylinder  is  zero  we  have 


T(X.t) 


[4Rap/n] 


n=<» 

53  C(-l)n/sin[(2n-l)irat/Ljsin[(2n~l)rX/L] 
n*=l  (5.15) 


[4RauL/ir  ] 


II— ~ ^ 

53  C C- 1 ) 1 / (2n-lj  ]sin[  (2n-  l)iratL]cos[  (2n-l) 
n=l 


X 

a/L] 

L/2 


(5.16) 


Using  the  fact  that  the  cosine  terms  vanish  when  X=L/2  and  the  definitions 

S = T/ur2  B stress  (5.17) 

P = M/(rrr2L)  (5.18) 
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1 


U = M/L 


(5.19) 


where  M is  the  mass  of  the  cylinder  and  r is  its  radius,  we  have  finally 

S(X,t)  = [4  Rap  li/nj  53  C C-D  /(2n-l)z]sin[(2n-l)Trat/L]cos[(2n-l) 

n=l 


nX/L] 


(5.20) 


On  the  basis  of  Equation  5.20  a number  of  conclusions  may  be  drawn. 
At  small  times,  i.e.  when 


sin[(2n-l)irat/L]  [(2n-l)irat/I,) , 


(5.21) 


Equation  5.20  reduces  to 


? n=“” 

S(X,t;  = [4pRa/ir  ] £ C(-l)n+  /(2n-l)](ira/L)tcos[(2n-l)TTX/L].  (5.22) 


If  we  use  the  series  identity. 


fl/4  =5^[(-l)n+1/(2n-l)]cosL(2n-l)TrX/L], 


Equation  5.22  reduces  to 


S(X,t)  = (pRa  t) 


(5.2i) 


(5.24) 


Equation  5.24  is  an  expression  of  our  expectation  that  a uniform  velocity 
gradient  produces  a uniform  stretching  until  a relief  wave  can  propagate 
from  the  ends  of  the  cylinder.  The  time  required  for  the  uniform 
stretching  to  exceed  the  elastic  limit  of  copper  in  our  "standard  jet" 
is  estimated  from  Equation  5.24  as 


t - S/(pRa  ) 


'*■  0. 16xl0"^sec. 


(5.25) 


We  see  that  the  strain  in  a stretching  jet  is  sufficient  to  insure  plastic 
flow  from  jet  formation  to  break-up.  Plastic  flow  reduces  the  velocity  a 
at  which  signals  propagate  so  that  what  happens  at  one  point  in  the  jet 
docs  not  communicate  rapidly  to  neighboring  points.  Therefore  the  jet 
continues  to  stretch  in  a manner  determined  by  the  local  velocity 
gradient  or  strain  rate. 


_iJ 


..AaTrcffc  iSiimunt^a 


If  the  deformation  of  the  cylinder  were  totally  elastic,  we  would 
expect  maximum  strain  when  every  term  in  Equation  5.20  is  positive  or 


i.at/L  = it/2,  (5.26) 

or 

t = L/2a.  (5.27) 

This  is  the  time  for  a sound  wave  to  travel  from  the  end  to  the  center 
of  an  elastic  cylinder.  In  this  case  Equation  5.20  becomes 

? n=°°  5 

S(X,L/2a)  = [4pRaL/iT]  £ [ 1/ (2n-l)  ]cos[ (2n- 1) ttX/l]  (5.28) 

n=l 


We  now  use  the  fact  that  the  Fourier  series  for  [(L/2)-X]  in  the  range 
-L/2  <X<L/2  is  given  by 


L/2-X  ■ (4L/t>2)  £ [ 1 / (2n- 1)  2]cos[  (2n- 1)  irX/L]  , (5.29) 

n=l 

so  that  Equation  5.28  becomes 

S(X,L/2a)  = pRa[L/2  - X]  (5.30) 


This  agrees  with  our  boundary  condition  that  the  stress  vanishes  at  the 
ends  and  our  expectation  that  the  stress  is  a maximum  at  the  center  of 
the  cylinder. 

The  preceding  discussion  of  the  jet  as  an  elastic  right  circular 
cylinder  illustrates  a number  of  points  which  are  expected  but  which  are 
perhaps  more  believable  if  demonstrated  by  a standard  mechanical  analysis. 
The  tension  along  the  jet  is  uniform  (Equation  5.25)  until  the  elastic 
limit  is  exceeded,  and  is  expected  to  remain  so  until  the  jet  begins  to 
break-up,  if  perturbations  are  expected  to  propagate  poorly  in  the 
plastically  stretching  jet.  With  a given  velocity  gradient,  the  maximum 
stress  available  (Equation  5.30)  is  proportional  to  the  product  of 
length,  velocity  gradient,  wave  velqcity,  and  density.  This  agrees  with 
our  earlier  calculations,  based  on  kinetic  energy  considerations,  that 
longer  segments  will  stretch  more  because  they  have  more  energy  per 
unit  mass  (Equation  4.14).  Since  the  tension  is  greatest  at  the  middle 
of  the  jet  (Equation  5.30)  a relieved  segment  (one  that  has  lost  its 
velocity  gradient)  will  be  stretched  more  in  the  middle  and  have  a 
tendency  to  be  bigger  at  the  ends,  where  stretching  is  least. 
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If  we  assign  a maximum  stress  level  s,  Equation  5.30  permits  a 
calculation  of  the  distance  d from  the  end  at  which  this  stress  in  exceeded 
and  plastic  deformation  begins: 

d s (L/2-X)  = S/pRa 

= 2xl08/(8.93«103*104*3750)  = 0.6mm  (5.31) 

This  means  that  segments  above  about  1mm  will  use  almost  all  of  their 
kinetic  energy,  relative  to  the  center  of  mass,  in  causing  plastic 
deformation. 


As  a matter  of  interest,  whether  or  not  it  has  physical  reality  in 
the  case  of  the  actual  jet,  the  relief  wave  starting  from  the  free  end 
of  the  jet  should  cause  the  jet  to  stretch  less  in  the  region  near  the 
end  and  hence  insure  that  the  tip  particle  will  have  the  largest 
diameter  and  will  have  a greater  length  than  average. 


If  we  imagine  a jet  to  start  as  a right  circular  cylinder  and  to 
-4 

break-up  in  10  seconds  with  open  spaces  of  one  cm  roughly  equal  to  the 
length  of  the  particles,  we  can  estimate  crudely  the  magnitude  of  the 
forces  involved.  The  jet  has  roughly  0.6  gms/cm  mass,  and  a mass  of  0.3  gm 
must  move  a distance  of  0.25  cm  on  the  average  to  vacate  a space  of  0.5  cm. 
Using  Newtons  second  Law,  we  derive  the  formula, 

F = 2MS/t2  (5.32) 

where  F is  the  average  applied  force  which  must  be  exerted  to  move  a 
mass  m a distance  X in  time  t.  Using  the  above  numbers. 


F = 2*0. 3xl0"3x0. 5x10" 2/(10'4)2  = 300n 


(5.33) 


The  surface  tension  v of  molten  copper  is  1.103  n/m,  so  that  its 
contribution  to  the  tension  T in  a jet  of  3 mm  diameter  D is  computed 
to  be 

T = irDv  * O.Oln.  . (5.34) 

4 

This  means  that  the  surface  tension  is  too  small  by  a factor  of  3x10 
to  produce  the  observed  motion,  and  that  the  jet  is  unlikely  to  have 
broken  up  because  it  was  a liquid  cylinder  unstable  under  forces  of 
surface  tension.  On  the  other  hand,  the  strength  of  copper  is  roughly 

2x10  n/m  , so  that  the  tension  in  the  jet  if  it  is  a cold  solid  is 

1.4xl03n.  A 20%  variation  in  this  tension  would  be  required  to  produce 

-4 

the  specified  displacements  in  10  secs,  so  that  even  a sizeable 
instability  might  not  be  sufficient  to  produce  the  necessary  action 
unless  it  grew  rapidly  from  the  beginning. 
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6. 


A SIMPLE  STOCHASTIC  MODEL  FOR  JET  BREAK-UP 


With  the  preceding  theoretical  framework  it  is  now  possible  to 
discuss  jet  break-up  in  a more  quantitative  fashion.  It  will  be 
assumed  that  the  copper  jet  stretches  as  a solid  until  enough  work  is 
done  on  it  to  cause  it  to  break  into  pieces.  The  distribution  in  length 
of  the  pieces  will  be  determined  by  the  energy  available  to  stretch  a 
piece  further  after  it  is  formed,  the  distribution  of  weak  points  or 
inhomogenities  in  the  jet,  and  the  properties  of  the  material  in  general. 

The  first  part  of  the  determination  is  given  by  Equation  4.14, 


E/[M(t/L)]  = (l/24)(Us/L)  V,  (6.1), 

which  relates  the  specific  kinetic  energy  relative  to  the  center  of  mass 
of  a segment  to  the  velocity  gradient  (Us/L)  and  length  £.  We  are 

justified  in  keeping  Equation  6.1  while  the  jet  stretches  and  begins 
to  break-up  because  we  have  shown  following  Equation  4.11  that  only  about 
7%  of  the  jet  energy  is  consumed  in  stretching  and  have  argued  that  the 
jet  stretches  uniformly  with  little  communication  up  and  down  the  jet 
because  of  the  rapid  rate  of  plastic  deformation. 

The  second  part  of  the  determination,  namely  the  distribution  of  weak 
points  in  a jet,  is  something  easy  to  believe  but  difficult  to  demonstrate, 
other  than  by  noting  that  the  jet  does  break-up  into  a somewhat  irregular 
pattern. 

Referring  to  Figure  3,  suppose  that  the  probability,  that  any 
segment  of  length  AX  will  break  during  the  time  interval  between  t 
and  t+At,  is  given  by  the  expression, 

P(t)  AX  AE  = B(t)  AX  AE,  (6.2) 

where  AE  is  the  work  done  on  the  segment  by  stretching.  The  assumptions 
contained  in  this  equation  are  (1)  that  the  probability  of  a weak  point 
being  in  the  segment  is  proportional  to  its  length  and  (2)  that  the 
probability  of  its  breaking  is  proportional  to  the  work  done  on  it  by 
stretching. 


Since  work  proceeds  as  a function  of  time  as  the  jet  stretches* Equation 
6.2  can  be  written  in  the  more  convenient  form, 

p(t)  AX  At  = b(t)  AX  At,  (6.3) 


where 


t <t<t  . 
o m 


w III 

Here  tQ  is  the  time  when  break-up  starts  and  t 


m 


(6.4) 

is  the  time  when  it  ends. 
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Figure  3.  Calculation  of  the  probability  that  n segments  of  total 
length  L will  hold  together.  Breaks  arc  postulated  for 
segments  i and  n,  and  the  probability  that  no  other  breaks 
occur  before  a relief  wave  arrives  is  calculated. 


The  inclusion  of  the  limits  on  the  time  t reflects  the  experimental 
observation  that,  in  stretching  a number  of  like  samples  of  copper,  the 
elongations  at  which  breaks  occur  are  bunched  between  upper  and  lower 
limits  which  are  rather  loosely  spaced.  Between  these  limits,  the  work 
done  is  roughly  proportional  to  the  elapsed  time  At.  The  parameter  b(t) 
is  characterized  by  the  material  and  its  condition  at  time  t. 

Suppose  that  the  time  interval  between  t and  t.  is  divided  into 

01 

M equal  intervals  of  length  At  and  that  the  jet  is  also  divided  into 

imaginary  segments.  The  probability  Pik  that  the  ith  segment  of  length 

(AX)^  will  survive  the  Kth  time  interval  (At).,  without  a break  is  given 
from  Equation  6.3  as  * 

Pik  - 1 - b(t)  (AX).  (At)K,  (6.5) 


i.e.  the  survival  probability  p„  is  unity  minus  the  probability  of  a 
break,  given  by  Equation  6.3, 


The  total  survival  probability  of.  the  segment  is  the  product  of  all 
the  individual  survival  probabilities,  i.e.. 


pi(t) 


Ll-b(t)(AX).(At)K] 


(6.6) 


Taking  the  logarithm  of  both  sides  yields  the  result. 


Since 


Hn  p.(t) 


in[l-b(t)  (At).  (At)K] 


iim  in  (1-y)  = -y 
y-K) 


(6.7) 


(6.8) 


and  the  products  b(t) (AX) . (At)„  will  become  arbitrarily  small  as  (AX). 

IK  1 

and  (At)^  approach  zero,  it  is  a good  approximation  to  write 


K=M 


*n  p. (t)  = - £ b(t)(AX).(At)„. 
1 K=1  IK 


(6  9) 


The  identification, 


K»M  t . 

£ b(t) (At) K = 

K=1  K . 


/: 


(6.10) 


b(t)dt, 


permits  Equation  6.9  to  be  written  as 


*n  pt(t)  - -(AX)i 


b(t)dt 


(6.11) 


Equation  6.11  is  a statement  of  the  probability  that  the  ith  segment 
of  length  (AX)^  will  survive  if  a relief  wave  arrives  at  time  t..  It 

is  assumed  that  the  process  of  jet  break-up  is  as  shown  in  Figure  4 
and  that  a segment  is  subject  to  breakage  until  it  is  relieved  by  a 
relief  wave  propagating  from  a nearby  break. 


In  a like  manner  it  can  be  stated  that,  if  a rod  of  length  L is 
to  survive  until  time  t^,  each  sub-element  (AX),  must  be  relieved 

before  that  time.  Suppose  that  the  rod  is  divided  into  n segments  of 
equal  length,  as  shown  in  Figure  3,  and  that  the  rod  is  formed  by  a 
break  on  the  left  at  time  tj  and  a break  on  the  right  at  time  t2. 

Suppose  that  a relief  wave  travels  from  the  left  break  with  a speed 
U (X)  and  arrives  at  the  ith  segment  at  time  t^  The  survival  probability 

of  the  ith  segment  is  given  from  Equation  6. 11  as 


Pi(ti)  = EXP[-(AX)i  b(t)dt 


•f 


(6.12) 


where 


,-vf 


dX/il  (X)  . 


(6.13) 


The  distance  X^  is  measured  from  the  left  end  of  the  rod  to  the  ith 

segment,  and  U(X)  is  the  speed  at  which  the  relief  wave  travels  from 
the  left  break.  The  integral  in  Equation  6,13  is  just  the  time 
required  for  the  relief  wave  to  travel  from  the  left  break  to  the  ith 
segment.  In  general,  the  relief  waves  from  the  right  and  left  will 
meet  at  some  time  t and  at  some  point  X ^ such  that 


^12 

bu  " h * 

J 0 


•r 


dX/U(X)  - t2  +|  dXVU(X'), 


(6.14) 


X12  5 L " X12 


(6.15) 
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where 


b 

I 


v//jy///A 


4 3 

nzm 


M i 


Figure  4.  Brea);  up  of  & jet.  Break  points  are  numbered.  Relieved 
areas  are  shown  by  cross  hatch.  Fach  segment  continues 
to  stretch  until  it  is  either  broken  or  relieved  by  a 
relief  wave  traveling  from  a free  end. 


i 
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From  the  above  considerations  and  reference  to  Figure  3,  the 
probability  thdt  a break  will  occur  on  the  left  between  the  times  t^ 

and  t^  + (At] 1 and  on  the  right  between  the  times  t.,  and  t2  + (At)2 

and  that  the  intervening  segments  will  survive  until  reached  by  a 
relief  wave,  after  which  they  are  sure  to  survive,  can  be  written  as 


P(L,t1,t2)(At)1(At)2  = 


K=n 


TTPk  b(tl)(b(t2)(At)l(At)2  (6>16) 


K=1 


In  the  above  expression,  b(tjHAt)j  is  the  probability  per  unit  length  that 

the  first  segment  will  break  between  the  time  t1  and  t^  + (At)^  and 

b (t2) (At) 2 is  the  probability  per  unit  length  that  the  last  segment  will 

break  between  the  times  tn  and  + (At)2>  For  F.quation  6.16  to  be  a 

true  probability  equation,  it  should  be  multiplied  on  both  sides  by 
(AX) 1 (AX) 2 to  include  the  lengths  of  the  segments  where  the  breaks 

occur.  However,  because  only  a relative  probability  is  sought  and 
because  the  location  of  the  breaks  is  immaterial  to  the  extent  that 
any  other  two  breaks  separated  by  the  same  distance  would  yield  the 
same  results,  the  factor  (AX)1(AX)2  is  ignored  to  save  writing  and 

space.  The  first  product, 

Y | P j > ' (6.17) 

M 

is  the  probability  that  the  left  part  of  the  rod  of  length  L,  defined 

by 

(XX<X12,  (6.18) 

will  survive,  and  the  second  product  is  the  probability  that  the 
right  part  will  survive.  Total  survival  probability  is  the  product 
of  the  survival  probabilities  of  the  individual  part3. 

For  convenience  of  analysis,  the  left  hand  side  defined  by  Equation 
6.18  is  divided  into  M equal  segments  of  length  AX  and  the  right  hand 
side  is  divided  into  n equal  segments  y length  AX  so  that 

(M+n)AX  =*  L (6.19) 
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From  Equations  6.12  and  6.13, 


1 

H 

pj s 

EXPL-AX  I 

b(t)dt]; 

t 

(6.20) 

o 

/k 

PK  = 

EXPL-AX  I 

b(t)dt]; 

^0 

(6.21) 

1 

rxj 

t.  = 
3 

l 

dy/U(y)  = tj  + f(X.); 

(6.22) 

rL-: 

h - 

*2  +J 

1. 

dyVU(y')  = t2  ♦ f(L-XK). 

(6.23) 

By  way  of  explanation,  Equation  6.20  states  the  survival  probability 
of  the  jth  segment  relieved  from  the  left  by  a relief  wave  starting  some 
;ime  between  t^  and  tj  + (At^  and  arriving  at  time  t j . Equation  6.21 

computes  t^  in  terms  of  the  initiation  of  the  relief  wave  at  time  tj 

and  its  travel  time  to  the  jth  segment  located  at  X„  Combining 

Equations  6.16-6.23  and  taking  the  logarithm  gives  the  result. 


in  P(L,t1(t2)  * tn[b(t1)b(t2)]  - AX 


2 (tj 
3=1  Jo 


K=n  f K 

b(t)dt  - AX  J2  I b(t)dt.  (6.24) 
K=1  J 0 


In  the  limit  as  AX  becomes  arbitrarily  small 

»Xi„  -t(tj,X) 

b(f.)dt  =|  dX  | b(t)dt,  etc, 

o o ' o 


A X 


i*M  f* j fhz  rt(tr 

22  b(t)dt  * I dX  I 
j*l  •'t  J t Jt 


(6.25) 


so  that  Equation  6.24  can  finally  be  written  as 


Jin  P(L,tj,t2)  * JUi[b(t 


(*12 

.)b(t2)]-Jo 


f « f 

Jo  J t 


(t2,X) 


b(t)dt 


/t  (t  j , X) 

b (t)dt 


(6.26) 
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The  integration  for  the  left  side  extends  from  time  tQ,  when  breaks 
can  begin,  to  the  time  t(t^,X)  when  the  relief  wave  arrives,  where 
from  Equation  6.22 

rx 

t(trX)  - ti  + dy/U(y)  . (6.27) 


Likewise, 


fL-X 

t(t2,X)  - t2  + dy VU(/)  . 


(6.28) 


The  integration  over  X,  for  example,  continues  from  the  left  end  to 
the  point  where  the  relief  waves  meet.  Equation  6.14  determines 

the  meeting  point.  Finally 


P(L) 


■//**.  dt2  P(L,t1t2).  (6.29) 


To  exercise  this  theoretical  structure,  it  is  worthwhile  to  make  the 
simplest  possible  assumptions,  namely  that 


b(t)  = b, 


(6.30) 


where  b is  a constant  and  that  relief  waves  travel  with  a constant  speed 
a,  i.e. 


U(X)  = a. 


(6.31) 


Equation  6.14  then  becomes 


/X12  fL'X  12 

dX/a  = t-  + I 

0 2 Jo 


dX/a  (6.32) 


- t:  ♦ X12/a  = t2  + (L-X12)/a  (6.33) 


so  that 


X12  “ ta(t2  ' tl)  + L]/2> 


(6.34) 
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| From  Equations  6.27  and  6.28, 

t(t1,X)  = ^ + X/a;  (6.35) 

t(t2,X)  = t2  + (L-X)/a.  (6.36) 

With  these  definitions  the  integrals  of  Equation  6.26  can  be  written 
as 

rX12  fX12 

I dX  J b(t)dt  = b I [t.  + X/a  - t ]dX 

' 0 J t J0 

o 


= b C(tj  t0)X12  + X^/2a] 

(6.37) 

and 

fL"X12  |^(t2iX)  rL"X12 

1 dX  I b(t)dt  » b 1 [t  +X'/a-t  ]dX' 

j 0 Jt  J 0 

o 

= b [(t2-to)(L-XJ2)  + (L-X12)2/2a] 

(6.38) 

After 

these 

some  simplification  and  the  use  of  Equation  6.34,  the  sum 
integrals  can  be  written  as 

of 

(ba/4)  [2T(X+y)  - (y-X)2  + T2] 

(6.39) 

where 

X 3 trtQ; 

(6.40) 

y = t2-to; 

(6.41) 

T s L/a. 

(6.42) 

Equation  6.29  can  now  be  written  in  the  form. 

P(L)  *=  b 2J  JdX  dy  EXP (ba/4)  [2T(X+y)- (y-X) 2+T]| 

(6.43) 
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A number  of  conditions  must  be  met  for  the  limits  of  integration 
indicated  in  Equation  6.43.  Since  the  break-up  occurs  between  t and 

t 

m 

W*.  ’ "‘“W*  '6'44> 

t0*'t  2*'tm  • (6-4S) 

Because  all  the  rod  of  length  L must  be  relieved  before  time  t , 

t<t.  (6.46) 

u m 

From  Equation  6.33  it  can  bo  shown  that 

tu  = [L/a+t1+t2]/2,  (6.47) 

so  that,  from  Equations  6.46  and  6.47 

L/a  + t,  + t0<2t  , T + X + y<2T.  (6.48) 

12m 

Finally,  the  breaks  at  times  tj  and  t2  must  be  such  that  a relief  wave 
from  one  does  not  arrive  in  time  to  prevent  the  other,  i.e. 

|t1-t1|<L/a  , |X-y|«?T.  (6.49) 

These  conditions  have  been  plotted  in  Figures  5 and  6 for  t^  and  t^. 

Since  there  is  symmetry  with  respect  to  which  break  occurs  first,  the 
arbitrary  condition 

t2<tj  (6. SO) 

was  added  to  reduce  the  range  of  integration. 
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Figure  6.  Graphical  representation  of  the  conditions  of  Equations  6,44*6.50 
for  l<a(tm-to).  The  cross  hatched  area  meets  all  the  conditions. 
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Unfortunately,  the  integral  of  Equation  6.43  cannot  be  expressed 
in  simple  terms.  However,  an  upper  limit  can  be  assigned  by  noting 
that  the  exponential  is  always  less  than  unity  so  that 

P(L)<b2A(L),  (6.51) 

where  A(L)  is  just  the  area  over  which  the  integration  is  to  be  performed. 
By  inspection  of  Figure  5 it  is  seen  that  the  area  of  integration  for 
I>a(tm"t0)  is  one  fourth  of  the  square  with  edge  of  length  [2 Ctm-tQ)  - 

L/a] . Therefore 

A(L)  - (l/4)[2(tm-t0)-L/a]2.  (6.52) 

In  Figure  6,  the  area  of  integration  is  seen  to  be  a trapezoid  with 
altitude  /?L/2a  and  bases /?(tm- L/2a-tQ)  and  /?(tm-L/a-to) . Therefore 


A(L)  ■ (l/2)(B^B2)h  = (t^t^L/a  -3L2/4a2)  (6,53) 

Combining  Equations  6.52  and  6.53  gives  the  result, 

ACL)  - (tm-to)CL/a)-3L2/4a2  , (KL-Laft^ t())  (6.54) 

■ U/4)[2(t,-t0)-L/a]2  , . 

It  is  seen  that  A(L)  has  a maximum  at 

L « 2a(tm-to)/3  (6.55) 

and  that  A(L)  meets  the  requirements  that 


A(0)  = A[2a(tw-tc)]  = 0,  (6.56) 

as  well  as  the  requirement  for  continuity  at  A[a(t  -tQ)]. 

The  maximum  for  P(L)  should  lie  to  the  left  of  the  maximum 
for  bZA(L)  because  of  the  increasing  influence  of  the  exponential  for 
larger  values  of  L. 


c 

An  experimental  plot  of  P(L)  is  shown  in  Figure  7.  A graph  of  A 
versus  L is  shown  in  Figure  8.  From  Equation  6.56,  if  the  maximum 
permissible  length  is  determined  from  Figure  7 to  be  3.8  cm  and  if  the 
relief  velocity  is  assumed  to  be  3750  m/sec, 


2a(t  -t  ) * L (6.57) 

m o max 

t -t  = 3.8* 10-2/ (2X3 750)  = 5X10~6  sec.  (6.58) 
m o ' 

The  ratio  of  maximum  length  tp  most  probable  length  is  roughly  3 from 
the  experimental  data  and  is  also  3 from  Equations  6.54  and  6.55. 

The  shape, of  the  theoretical  curve  is  not  very  good. for  small 
values  of  L,  possibly  because  of  the  variation  of  a with  distance 
from  the  break.  The  relief  speed  a is  probably  near  375C1  m/sec^at 
the  break  but  less  as  the  energy  dissipation  required  of  the  relief 
process  is  increased,  according  to  Equation  4.14, 

The  limits  of  integration  for  Equation  6.29  can  be  determined  from 
Figures  5 and  6 by  inspection.  From  Figure  5,  for  I>a(tm-to), 


t.<t.<2t  -t,-L/a;  (6.59). 

12  m i 

t <t,<t  -L/2a.  (6.60) 

o 1 m 

From  Figure  6,  for  0<L<a(tn)'  to)  and  t^t^t^-l/a, 

t1<t2<L/a  + tjj  (6.61) 

t <t,<t  -L/a.  (6.62) 

o 1 m 

From  Figure  6,  for  0<L<a(tffl-to)  and  tm-L/a<t j<tn~L/2a, 

t.<t_<2t  -t.-L/aj  (6.63) 

12  mi 

t -l/a*St.<t  -L/2a.  (6.64) 

m t m 


^Robert  L.  Jameson  of  the  Ballistic  Research  Laboratory,  private 
communication,  Nov  1977 
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These  limits  can  be  changed  to  suit  Equation  6.43  by  the  substitutions 
of  Equations  6.40-6.42. 


For 

T»tHt  -t  , 
m o 

(6.65) 

1 

r 

X<y<2t-X-T, 

(6.66) 

0<X<t-T/2. 

(6.67) 

For  Kt. 

x<y<r+x, 

(6.68) 

0<X<t-T; 

(6.69) 

X<y<2t-X-T; 

(6.70) 

. t-KX<t-T/2. 

(6.71) 

With 

these  limits,  Equation  6.43  may  be  written  more 

/•t-T  /T+X  /.t-T/2  «2t-X-T 

explicitly  as 

P(T) 

- | dX  1 F(X,y)dy  ♦ f dX  | 

J 0 j X j t-T  j X 

P(X,y)dy  (6.72) 

when 

0*2Kt  and  as 

/t-T/2  -2t-X-T 

dX  I F(X,y)dy  (6.73) 

0 J X 

when  t<K2t,  where 


F(X,y)  = EXP  |-ba[T(X+y)/2  + T2/4  - (X-y)%]|*  (6.74) 


7.  FURTHER  CONSIDERATION  OF  JET  BREAK-UP 


It  may  seem  reasonable  from  a physical  viewpoint  to  remove  the 
time  restrictions  on  the  onset  and  completion  of  jet  break-up.  Experience 
with  pulling  copper  seems  to  justify  the  assignment  of  a narrow  range 
between  the  strain  necessary  to  start  breakage  in  a collection  of  samples 
and  the  strain  sufficient  to  insure  breakage  in  every  sample.  Experimentally, 
there  arc  no  samples  which  do  not  break,  and  it  is  not  reasonable  to 
expect  that  a segment  of  a jet  will  be  arbitrarily  long  after  jet 
break-up  and  relief  of  the  individual  segments. 

The  process  by  which  a segment,  with  a velocity  gradient  and  two 
free  ends,  relieves  itself  is  not  well  understood.  If  it  is  assumed, 
as  shown  in  Figure  9,  that  the  fragment  of  length  L is  relieved  from 
the  left  end  to  the  point  X,  which  has  a speed  RX,  the  kinetic  energy 
Er  for  the  relieved  part  can  be  written  as 

Er  = (M/2L)(L/2-X)R2X2  (7.1) 

The  kinetic  energy  Eu  of  the  unrelieved  part  is  expressed  as 


R2X2dX  = (M/6L)R2X3  (7.2) 


The  total  energy  E of  the  left  half  of  the  segment  is  therefore  given 
by  the  equation, 

E =■  Ed+E  « MR2X2/4-MR2X5/3L.  (7.3) 

K U 

It  will  be  assumed  that  this  energy  is  removed  by  stretching  of  the 
unrelieved  portion  of  the  segment.  The  rate  P at  which  work  is  done 
in  stretching  the  unrelieved  portion  is 


Eu  = (M/2L) 


/ 
J 0 


P « TRX  (7.4) 

where  T is  the  tension  at  the  interface  between  the  relieved  and 
unrelieved  regions  and  RX  is  the  velocity  of  the  material  at  the  inter- 
face. In  an  instant  of  time  dt,  the  work  done  dE  can  be  written  as 

dE  = Pdt  = TRXdt . (7.5) 
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The  energy  loss  by  the  segment  is  given  by 


(:&E/3X)dX  = (MR2X/2-MR“X2/L)dX  « dE  (7.6) 

where  E is  given  by  Equation  7.3  and  dX  is  the  amount  the  interface 
moves  in  time  dt.  Equating  tho  results  of  Equations  7.5  and  7.6 
gives  the  equation, 

= (TRX)/'(MK2X/2-MR2X2/L)  (7.7) 

= SAL/ (MRL/2-MRX)  =■  (S/pR)  (L/2-X)  (7.8) 

where  the  tension  T has  been  written  in  terms  of  the  stress  S and 
the  cross  section  C of  the  segment  as  SC  and  the  density  p has  been 
defined  as  M/LC. 

If  t.he  defination, 

X'  * (L/2)  - X (7.9) 

is  made.  Equation  7.S  becomes  ' i 

dX/dt  - S/pRX  '=1/ (AX')  (7 . ill)  1 

1 

This  is  the  speed  which  the  interface  (relief  wave)  must  move  to  < 

consume  the  available  energy  if  the  segment  is  relieved  as  hypothesized  1 

in  Figure  9,  where  X"  is  the  distance  from  the  left  of  the  segment.  j 

Following  the  pattern  of  Section  6,  it  is  seen  that  the  time  for  j 

the  relief  wave  to  arrive  at  a point  X from  the  end  is  given  by 

dX/lJCX)  = t.  + AXdX  - t.  + AX4/ 2 

0 Jo  . 

Likewise, 

t(t2,X')  - t2  + AX'2/2. 

The  waves  meet  at  a time  ty  and  place  X12  such  that: 

tx  + AX22/2  =>  t2  ♦ A(L-Xj2)2/2  « t 


t (t  j . X'j  * t,  * 


1 


(7.11) 


(7.12) 


(7.13) 
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so  that 


X12  " Ct2-tx)/(AL)  + L/2 

and  n n a 

tu  = (t2+t2)/2  + (t2-t1)V(2ALZ)  ♦ AL  /8  , (7.14) 

Proceeding  as  before  the  integrals  of  Equation  6.26  can  be  written  as 


fX12 

>.t(trX) 

(l~h2 

/•t(t2,X) 

| dX 

I dt  + b 1 

1 dX' 

1 dt,  (7.15) 

^0 

' 0 J 

' t 

or  more  explicitly  as 

1 = b{CWX12  + Ca/6)X12  ♦ (t2-tQ)(L-X12)  + (A/6)(L-X12)2|  (7.16) 

Thi s reduces  to  the  expression, 

I * b[(X+y)L/2  t AL3/24  - (y-X)2/2AJ.]  (7.17) 

where  as  before 


XBVto 

(7.18) 

y = t2"to 

(7.19) 

1 B V'c  • 

(7.20) 

The  integration  of  EXP(-I)  over  the  variables  t2  and  t2 
several  restrictions  of  the  same  kind  as  before,  namely, 

is  subject 

(7.21) 

r»  ■■  • - 1 . • 

(7.22) 

|t2-tj|^L2/2,  |y-X|<AL2/2ST/2, 

(7.23) 
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(7.24 


tm>tu  - (tj*t2)/2  ♦ (t2-t1)2/(2AL2)  + AL2/8  , 
2t>(X+y)  + (X-y) 2/T  + T/4. 


Using  the  definition 
Equation  7.17  can  be  written  as 


I - (bLT/2)  [ (X+y) /T  + l/12-(y-X)2/T2] 


(7.25) 

(7.26) 

(7.27) 


Equation  6.39  can  be  written  for  the  constant  relief  velocity  as 

I - (baT2/4)  [2(X+y)/T  ♦ 1 - (y-XZ/T2]  (7.28) 

These  functions  (Equations  7.27  and  7.28)  are  somewhat  similar  in 
appearance  but  in  one,  Equation  7.27,  T is  defined  as  AL2  whereas 
in  Equation  7.28  it  is  defined  as  L/a.  This  means  that  longer  segments 
are  discriminated  against  when  the  relief  wave  speed  varies  as  1/AX. 
Small  segments  are  also  discrimminated  against  because  the  relief  is 
very  rapid  near  a break  so  that  the  next  break  is  more  likely  to  be 
distant.  Thus  the  1/AX  assumption  for  the  relief  wave  speed  should 
result  in  a more  narrow  distribution  than  the  assumption  of  a constant 
relief  wave  speed. 

To  illustrate  the  above  points  more  specifically,  from  conditions 
described  by  Equations  6.47  and  6.48  and  examination  of  Figure  6,  it 
is  noted  that,  as  L approaches  zero,  the  trapezoid,  which  represents 
the  area  over  which  the  and  t^  integrations  take  place,  approaches 

a height  of  l/(/5a)  and  a length  /2(t-L/a) . Therefore  as  L approaches 
zero 

Area  ■ Lt/a  (7.29) 

This  behavior  at  small  L is  exhibited  by  the  curve  of  Figure  8 near  the 
origin. 

On  the  other  hand,  for  the  1/AX  hypothesis  for  the  relief  wave  velocity, 
the  condition  of  Equation  7.23  indicates  that  the  area  of  integration, 
corresponding  to  that  of  Figure  6,  has  a width  of  AL2/2.  The  length 
of  the. area,  from  Equation  7.24,  in  the  limit  as  tj  and  t^  are  equal 

and  approach  tg,  is  /J  (t-AL2/B).  In  the  limit  as  L goes  to  zero,, 
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Area 


(7.30) 


AL2//f. 


This  behavior  is  characteristic  of  the  experimental  curve  of  Figure  7 
near  the  origin.  The  hypothesis  of  a 1/AU  relief  velocity  appears 
to  yield  a theoretical  curve  which  is  reasonable  in  shape  near  the 
origin  and  at  large  values  of  L,  but  the  exact  form  of  the  curve  is 
not  known  without  evaluating  the  integral  of  EXP(-l),  given  by  Equation 
7.16  for  various  combinations  of  the  parameters  b.  A,  and  t. 


8.  CONCLUSIONS 

A simple  physical  model  has  been  developed  to  explain  jet  break-up. 

It  gives  a fairly  good  description  of  the  upper  limit  for  particle 
length  after  break-up  ai'id  predicts  particle  distribution  with  respect 
to  length  in  a relative  way.  The  model,  unfortunately,  makes  use  of 
material  parameters  which  are  not  readily  observable,  such  as  the 
strength  of  a stretching  jet,  the  limit  to  which  the  material  can  be 
worked  by  stretching  without  breaking,  and  the  speed  of  a relief  wave 
originating  at  the  free  end  of  a segment.  It  is  probable  that  any  model 
would  have  the  same  problem,  because  these  parameters  are  fairly 
fundamental. 

Some  obvious  suggestions  for  improving  jets  arise  from  the  study. 

If  the  size  of  particles  has  an  upper  bound  associated  with  the 

development  of  breaks  at  propitious  times  and  places,  why  not  design  j 

the  jet  to  be  weak  at  intervals  whose  distance  is  slightly  less  than  •! 

the  maximum  permissible  length?  The  jet  could  then  consist  of  particles  > 

all  of  which  are  near  maximum  permissible  length.  ] 

J 

The  time  between  onset  and  completion  of  jet  break-up  is  estimated  ; 

to  be  SxlO  ^ seconds,  and  is  inversely  proportional  to  the  speed  of  the 
relief  wave  propagating  from  a break.  This  corresponds  to  a strain  of  1 

roughly  5%  between  onset  and  completion  of  break-up  for  the  jet  under  , 

study.  j 

If  the  jet  must  break  rather  than  stretch  indefinitely,  a more  j 

uniform  jet  will  break  into  smaller  pieces  than  a less  uniform  jet.  j 
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